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ABSTRACT 



Blue primordial power spectra with a spectral index ?T- > 1 can lead to a significant 
production of primordial black holes in the very early Universe. The evaporation 
of these objects leads to a number of observational consequences and a model 
independent upper limit of n ~ 1.4. In some cases this limit is strengthened to 
^1 77. = 1.3. Such limits may be employed to define the boundary to the region of 

parameter space consistent with generalized inflationary predictions. [To appear 
in Proceedings of the CASE WESTERN CMB WORKSHOP, April 22-24 1994. 
• Figures available on request from J.H.Gilbert@qmw.ac.uk] 
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G^ ■ 1. Introduction: Towards an Observational Test of Inflation 

^: 

Ph' The positive detection of anisotropic structure in the temperature distribution 

Q ■ of the Cosmic Microwave Background (CMB) radiation has opened up the possibil- 

ity that the predictions of the inflationary scenario may be testable within the near 
^ \ futurea'B. During inflation the scale factor grows exponentially, whilst the Hubble 

radius H^^ ~ lO^^^cm remains almost constant. Consequently the physical wave- 
length of a quantum fluctuation in the scalar or graviton fleld soon exceeds and 
its amplitude then becomes 'frozen'. Once inflation has ended, however, in- 
\ creases faster than the scale factor, so the fluctuations eventually reenter the Hubble 

radius during the radiation- or matter-dominated eras. The fluctuations that exit 
around 60 e-foldings or so before reheating reenter with physical wavelengths in the 
astrophysically interesting range 1 Mpc - 10^ Mpc. Fluctuations in the graviton de- 
grees of freedom result in a stochastic background of primordial gravitational waves 
that have an amplitude at reentry given by (5gw ~ ^/"^pb where m-pi is the Planck 
mass. Fluctuations in the inflaton fleld provide the seeds for galaxy formation via 
gravitational instability and their amplitude at reentry is 5 ~ V'^''^/(mpi|y |). Both 
scalar and tensor fluctuations lead to CMB anisotropics. 

Since H decreases as the inflaton fleld rolls down the potential, the amplitude 
of the scalar and tensor fluctuations is scale-dependent. This variation is most 
conveniently parametrized in terms of the spectral indices, ut and n, deflned by 
(5gw(M) cx M-^'T/a and 5{M) oc M(l-")/^ where M is the mass scale associated 
with the Hubble radius at the epoch of reentry. In general these spectral indices 



are themselves functions of scale. However, scales relevant to large-scale structure 
and CMB experiments probe only 9 e-foldings or so of the inflationary expansion, 
and since the scalar field must be rolling slowly for inflation to occur in the first 
place, these scales correspond to a very small portion of the inflationary potential. 
Hence, it is reasonable to suppose that the spectral indices are indeed constant over 
the scales of interests. 

When comparing theory with experiment, it is conventional to expand the CMB 
temperature fluctuations on the sky in terms of spherical harmonics. The angular 
correlation function predicted from theory is then given by an average over all 
observer positions: 

{6T{e)6m) = y: -^QP,{cose), (1) 

1=2 

where the P/'s are Legendre polynomials and a given multipole / corresponds to an 
angular scale 9/1" ^ 60/1. The C^'s corresponding to the scalar and tensor 

(C; r) fluctuations are determined once the precise functional forms of the spectra 
have been specifled. In the limit that |nr| and |?2 — 1| are constant and small, the 
ratio of the / = 2 multipoles is uniquely determined by the tensor spectral index: 

R^^^ -7nT. (2) 

^2,S 

This consistency equation is a fairly generic prediction of inflation!. The quan- 
tities in this expression are measurable, at least in principle, and it forms the basis 
for an observational test of the scenario. If no reionization occurs, experiments 
on scales 9 > 2° [l < 30) provide a measure of the sum of the scalar and tensor 
contributions: Ci = Ci^s + C*i,T- Since the gravitational waves do not produce a 
measurable contribution to the CMB anisotropy for 9 < 2°, one might hope to de- 
termine C2,T and 6*2,5 separately from a combination of small and large angle CMB 
experiments. A test of inflation would then require a separate determination of the 
tensor spectral index. At present it seems that this requires a direct detection of 
the gravitational wave spectrum. Recent calculations! suggest that the maximum 
present-day contribution per octave of the gravitational waves is Qcwh"^ < 7x 10^^^. 
This is too weak to be detectable by the Laser Interferometer Gravity- Wave Obser- 
vatories, although the proposed beam-in-space experiment has a peak sensitivity of 
Qqw ~ 10~^^ at 10~^ Hz and may be able to detect such a background. 



2. Constraining the Scalar Spectral Index with PBHs 

The scalar spectral index is much easier to measure and we would like to have 
an expression equivalent to the consistency equation that relates the C2's to n. In 
general, the relationship between and n is model dependent, but for the special 
case of an exponential potential, we have n — 1 = ut- Hence, in the parameter 
space {R,n), this model corresponds to the line R = 7(1 — n). It is important to 
note, however, that other families of potentials will lead to different trajectories in 



this space. For example, any potential of the form V = Vo[l ± 27i\n — ^(p'^Jmpy] 
leads to a measurable deviation of n from unity, whilst predicting that R ~ O.H One 
could imagine taking all the current candidates for the inflationary potential and 
predicting their separate trajectories in this parameter spacer. A superposition of 
these paths would then deflne a target that represented a generalized prediction 
of inflation in some sense. The problem with this approach, however, would be in 
deciding which models should be included in the target. Whilst a given model may 
appear to be natural to one person, someone else may deem it unnatural. 

It seems to us that this subjective element in the testing procedure must be 
eliminated if such an approach is to be developed further. To accomplish this it 
is necessary to search for model-independent constraints on the parameters R and 
n. Recently an upper limit on the spectral index was derived from considering the 
formation and subsequent evaporation of primordial black holes (PBHs)i. This limit 
is independent of any gravitational wave contribution to the CMB anisotropy and 
therefore deflnes a boundary to the observational target of inflationary predictions 
in {R, n) space. 

The idea behind the argument is rather simple. During inflation the flrst scales 
to leave the Hubble radius are the last to come back in and this implies that the very 
last fluctuation to leave is the flrst to return. In the simplest case, the fluctuation 
on this scale will be spherically symmetric and Gaussian distributed with an rms 
amplitude given by 5(te), where te ~ is the time when inflation ends. In some 
regions of the post-inflationary Universe, the fluctuation will be sufficiently large 
that the collapse of a local region into a black hole will become inevitable. The 
higher the rms amplitude the higher the fraction of the Universe forming PBHs. 
The observational consequences of the evaporation of these black holes then leads 
to upper limits on the number that may form and hence on the magnitude of 5{te)- 
An upper limit on n is therefore derived by normalizing the power spectrum on the 
quadrupole scale, Mq ~ lO^^g, and assuming that the spectral index is constant. 

PBHs are never produced in sufficient numbers to be interesting li n < 1, but 
they could be if the spectrum is 'blue' with n > 1. The precise form of the con- 
straints depends crucially on how the Universe is reheated, however, and we now 
proceed to discuss the separate cases of efficient and inefficient reheating. 

2.1. Efficient Reheating 

When an overdense region with equation of state p = stops expanding, it 
must have a size greater than y/j times the horizon size in order to collapse against 
the pressure and this requires that 5(te) > 7- It follows that the probability of a 
region of mass M forming a PBH isS 

A(Af)«*(M)exp(-^). (3) 

Because of the exponential expansion PBHs that form before or during inflation 
have no observational consequences. Furthermore, if the reheating process is very 



efficient, tlie false vacuum energy is rapidly converted into relativistic particles with 
a reheating temperature Trh/^pi = (^e/^pi)~^^^- The mass of a PBH forming at 
this time is Mrh/'^pi ~ %h/^pi and once it has formed, a PBH of this mass will 
evaporate at a time tevap ~ (^RH/"^pi)^tpi. Eq. then implies that for a blue 
spectrum, Po{M) decreases exponentially for M > Mrh, so we may regard the PBH 
mass spectrum as effectively being a 5-function at Mrh- 

The constraints on /3o(M) in the range lO^^g < M < 10^''g were recently 
summarized^. In particular, PBHs with an initial mass ~ lO^^g would evaporate at 
the present epoch and may contribute appreciably to the observed gamma-ray and 
cosmic-ray spectra at 100 MeV. On the other hand, 10^°g PBHs have a lifetime ~ 
1 sec and, if produced in sufficient numbers, their evaporations would lead to the 
photodissociation of deuterium immediately after the nuclesynthesis era. PBHs of 
mass slightly below could alter the baryon-to-photon ratio just prior to nucle- 
osynthesis. In our papeiCl we consider the constraints on Po{M) below 10^°g. In this 
region there is a potentially stronger constraint on the spectral index if evaporating 
PBHs leave stable Planck mass relicsEl. Although the formation of such objects has 
not been proved conclusively, it would be surprising if quantum gravity effects did 
not become important once the PBH had evaporated down to the Planck mass and 
various arguments have been developed in the literature suggesting the formation 
of such objects is likelyQ. 

To derive the observational constraint from PBH relics one proceeds as follows: a 
lower limit on Mrh is derived by assuming that the observed quadrupole anisotropy 
is due entirely to gravitational waves. This implies^ that the expansion rate of the 
Universe during the last 60 e-foldings of inflation cannot exceed 3 x 10~^mpi and 
leads to an upper limit on the reheat temperature of ~ 10^^ GeV. This corresponds 
to a minimum mass ~ Ig. The observational constraint from the relics derives from 
the fact that they cannot have more than the critical density at the present epoch, 
i.e. Qrci < 1- The precise form of the constraint depends on whether the evaporating 
PBHs dominate the energy density of the Universe before they evaporate. Since 
the ratio of PBH density to radiation density increases as t^/^ the PBHs do not 
dominate at evaporation if Pq < (Mrh/'^pi)"^- If this condition is satisfied, the 
constraint that flre\ does not exceed unity becomes 

/5o(M) < 10-2^ ( — ) . (4) 

If, on the other hand, the PBHs dominate the density at evaporation, most of the 
background photons derive from the PBHs and the constraint becomes M > lO^g. 
In other words, only relics formed from PBHs smaller than lO^g can contribute 
significantly to the current energy density and above this critical mass the entropy 
constraint takes over. 

We have completed a detailed analysis and derived the constraints on (5(Mph) 
for all mass scales above Ig. These results are illustrated in Figure 1. It is seen 
that the strongest limit on the spectral index derives from the relic constraint and 



by combining Eqs. (|^) and (||) one finds that 




.-1/2 



5(Mrh) < 0.13 17 -log 



(5) 



If we normalize on the quadrupole scale, Mg, the rms amplitude on a smaller scale 
M is 6{M) = 6q{M/MqY^-''^/^, where 6q ^ 3.8 x 10"^ Hence, we conclude from 
Eq. (1^) that the limit on the spectral index is n < 1.4 for Mrh ~ Ig? corresponding 
to a reheat temperature ~ lO^^GeV, and n < 1.5 for Mrh ~ lO^g, corresponding to 
^RH ~ lO^^GeV. The constraints associated with higher masses (i.e. lower reheat 
temperatures) are calculated by a similar procedure and are summarized in Figure 
2. 



Fig. 1. The constraints on the rms amplitude of the scalar fluctuation spectrum immediately after 
inflation if the equation of state is radiatipn-like (7 = 1/3). The origin of the constraints above 
10^'^g is summarized by Carr and LidseyS. If PBHs do not leave behind stable relics after evapo- 
ration, the strongest upper bound on the spectral index is given by the dashed line which joins the 
COBE/DMR point and the deuterium constraint at lO^^g. This limit applies for reheat temper- 
atures ~ lO^GeV. If relics are formed, the limit is strengthened at higher reheat temperatures as 
indicated. 



2.2. Inefficient Reheating 



Thus far we have assumed that the reheating of the Universe to relativistic 
particles occurs on a timescale much less than H~^. However, when inflation ends 
by means of a second-order phase transition, the scalar field undergoes coherent 
oscillations in the potential minimum from the time ti ~ H'^ until it decays at 
a time ^2 ~ T"^: where F is its decay width. During this interval, the Universe 
is effectively dominated by a dust-like fluidS with 7 = 0, and Eq. does not 
apply. PBHs will still form in this case, but the fraction of the Universe going into 
PBHs is now determined by the probability that regions are sufficiently sphericallv 
symmetric to collapse within their Schwarzschild radius. This fraction is given byt3 

(3{M) ^ 2 X 10-^[6{M)]^^^^ (6) 

and the observational constraints on the probability of PBH formation are altered 
because they now have an extended mass spectrum. The range over which this 
spectrum applies is defined by Mi < M < Mmax? where Mi is the horizon mass im- 
mediately after inflation and Mmax is the mass of a configuration that just detaches 
itself from the universal expansion at t2- It is determined implicitly bylll 

M^ax = [S{M^..)f' (^) mpi. (7) 

Vtpi/ 



The constraints on 13{M) are related to the associated constraints on I3q{M) via the 
relation 

Vtpi/ \m-p\J 

This is a useful expression because it implies that the limits on PBH formation 
during the dust phase can be calculated directly from the constraints which apply 
if there is no dust phase. We first assume ti is fixed and vary the epoch t2 = at 
which the dust era ends. For a given value of t2, the mass range of PBHs forming 
during the dust phase goes from Mi to the mass given by Eq. (|^) and for each 
value of M the constraint on j3{M) is given by Eq. (P). The corresponding limit 
on 5{M) follows from Eq. @ and the limit on n is derived by normalizing on the 
quadrupole scale as before. The reheat temperature is Trh ~ (r^pi)^''^'^pi and an 
upper limit on Ftpi follows from the requirement that baryogenesis must proceed 
after reheating. It is generally accepted that the lowest temperature for which the 
observed baryon asymmetry may be generated is the electroweak scale, ~ lO^GeV, 
corresponding to Ftpi ~ 10"'^°. 

The new constraints on the spectral index are also shown in Figure 2. For a 
relatively short dust phase, only the relic limit will be altered, since PBHs above 
10^°g will not form during the dust phase. For lower reheat temperatures, however, 
more massive PBHs form and provide the strongest constraint. The deuterium 
constraint applies for 10~^^ > Ftpi > 10^^^ and the gamma-ray limit for 10"^'^ > 



rtpi > 10 From this figure we arrive at an upper limit of n = 1.4 if PBHs form 
relics and there is a dust-phase after inflation. 



Fig. 2. Illustrating the constraints on the spectral index arising from the overproduction of primor- 
dial black holes, the shaded area being excluded. The lower line applies if there is a dust phase 



tic particles. The ?2-independent upper and lower limits on the decay width arise from assuming the 
COBE/DMR detection is due entirely to gravitational waves (shaded line) and from requiring that 
baryogenesis can only proceed above the electmweak scale (dashed line) . The dotted horizontal line 
indicates the CMB distortion limit of Hu et alli^. For reheat temperatures above Trh ~ 10^'^GeV 
the most important constraint arises from the requirement that any Planck mass relics left over 
from the final stages of PBH evaporation should have less than the critical density at the present 
epoch. For lower reheat temperatures, more massive PBHs may form and the strongest constraints 
then arise from the photodissociation of deuterium by evaporating 10^'^g PBHs, from the observed 
gamma-ray background in the energy range 0.1 — 1 GeV, or from the distortions of the CMB. 

3. Conclusion 

One criticism of these limits is that they assume the spectral index is constant 
over the full range of scales. To answer this point we note that the general inflation- 




ary potential leading to spectra with constant spectral index > 1 is a combination 
of trigonometric functions with a Taylor expansion of the formQ 



V = Vo l + 27r(n-l)-V • (9) 



m 



PI J 

Any potential that leads to n > 1 will have a Taylor expansion of this form. It 
follows that, as the field rolls towards the minimum, the approximation of Eq. (^) 
to the general trigonometric potential becomes more accurate and so the variations 
in the spectral index become smaller. Consequently, one need only show that the 
spectral index is effectively constant over the scales corresponding to large-scale 
structure (53-60 e-foldings from the end of inflation). It is straightforward to show 
that this is the casei. 

We conclude that the formation of PBHs from quantum vacuum fluctuations 
immediately after inflation constrains the scalar spectral index to be less than 1.4. 
It is important to emphasize that, because our limit spans the large range of scales 
from ~ Ig to ~ lO^^g, (for comparison, all large-scale structure measurements span 
only 10 decades of scale), it is essentially independent of the errors that arise in the 
COBE normalization from possible tensor contributions, cosmic variance and the 
effect of the Doppler peak on the low multipole anisotropies. Moreover, since we 
have normalized on COBE scales, the limit is also independent of the precise form of 
dark matter and hence the bias parameter. In effect, once the COBE normalization 
is specified, the limit is independent of the cosmological model, although it does 
assume that = 1, as predicted by most inflationary models. Therefore this limit 
can be employed to deflne a boundary to the target of inflationary predictions in 
{R, n) space. 
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